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ABSTRACT 



We construct spectral decomposition of 3D Fokker - Planck differential operator in 
this paper. We use the decomposition to obtain solution of Cauchy problem - and espe- 
cially the fundamental solution. Then we use the decomposition to calculate macroscopic 
parameters of Fokker - Planck flow. 



1. Fokker - Planck 3D differential operator 

The Fokker - Planck equation in three dimensions is 



where 

n = n(t, x k , v k ) - density; 

t - time variable; 

x k - space coordinate; 

v k - velocity; 

a - damping coefficient; 

k - diffusion coefficient. 

(1) is evolutionary equation with differential operator 




(1) 




(2) 



where tj> = <p(x k , v k ). 



We call differential operator (2) - Fokker - Planck 3D differential operator. This operator is not self- 
adjoint. Its adjoint operator is: 



We discuss simple cyclic boundary conditions 

<j>{Xi + a h v k ) = <p(x h v k ), - oo < v < oo. 

where a k - period along the corresponding coordinate axe. 

2. Eigenfunctions and eigenvalues of Fokker - Planck 3D differential operator 

Eigenfunctions <j> of differential operator (2) must satisfy with differential equation 



L(fl = ^=£ -v* 



dx k 



+ a 



d(» Vfc) 



+ k 



dv k dv k 




= 0. 



and boundary conditions (4). 



We search solution of equation (5) in the form 



<P = X 1 (x 1 )X 2 (x 2 )X 3 (X3)y i (v 1 )V2(V2)V 3 (V3). 



or 




The separation of variables provides two equations 




Wr + aivtViY-CiVtVi 



where c, and fi t are separation constants. 



According to (7) 



The solution of (8) and (4) is 
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x k 

X hn = e «*; (11) 



m k 

where m k is integer and accordingly c k = 2m — . Now (9) takes the form 

"A 



m 



kV/' + a(VjVjY - 2m — '- VjVj = fijVj. (12) 



a j 



We solved ID equation (12) in our work [1]. Eigenvalues of (12) are 



( 2itm j x2 



z-yt lit , i 

jUj - -arij - k ^ -|. (13) 



Corresponding Eigenfunctions are 



2mm ; ( a r\ ( \~a~( Anim.;k\\ 

v ^ = exp( - ^ exp [- 2, *>y »(v 2* b + ^7 j) (u) 



where: 

- 1 < j <3 - coordinate number; 

- tij > and TO; - integers; 

- //„(#) - Hermite polynomial. 

Hermite polynomials are defined as 



(-1)" d" 

H n (z)="—— P ; (15) 

yO dl n 



where p = e z . 

Eigenfunctions according to (6) have the form 



"m l m 2 m 3 nin 2 n2 



H ( m i Vi \ f a 7 \ ( ra~( 4xim;kX\ 
= S «P ^ ^ " i) J«P [- 21 ^>,[V 2k[ Vj+ ^a~)) (16) 

Each eigenfunction is uniquely characterized by a set of 3 nonnegative integers rij and a set of 3 inte- 
nt 

Expression for eigenfunctions of adjoint differential operator is according to [1] 



: n «p [-™ -i + i> J *-(V 2i b + ^7 J) (17) 



Evidently, these eigenfunctions are orthogonal 



CL\ £?2 a l OO OO OO 

| dXi | dx 2 | | <*. 

— oo — oo — OO 



„„„„„ „ „ : „ , <?/ /v l ^2^-2'" ^""p,!p 2 !p3! | —J x (18) 



xa l a 2 a 3 £ m2 „ 2 £ m3 „ 3 8 Piqi 8 P2 q 2 S P3 q 3 exp 



a (\nkX 
2k 



( 



m \\ I m 2 ) ( m 3 



v y 



We can write following expansion of arbitrary function G(x ■ v.) 



+00 +00 +00 



G(Xj, V j) 2_j ^ ^ ^ ^ J ^m l m 2 m 3 p l p 2 p 3 ^tn l m 2 m 3 p l p 2 p 3 ^ 

mi=-oo m 2 =-oo m 3 =-oo /7j=0 /?2-0 P3=0 



(19) 



where coefficients A are equal to (see (18)): 



a X 



2 1 



exp 



2fc ["^"J 



mi m 2 m 3Pl p 2 p 3 2 Pt + P2+Pip l \p 2 \p J \{2zk) a x a 2 a 3 

&\ &2 a 3 OO OO OO 

dxx ^ dx 2 ^ dx 3 j" J JgC^v,) Wm 1 m 2 m 3PlP2 p s dv l dv 2 dv 3 . 



— oo — oo — OO 



x(20) 



3. Solution of Cauchy problem and fundamental solution 



Solution of Cauchy problem with initial condition n Q = G(x p Vj) is 



+oo +00 +00 oo oo oo 

n(t,Xj,Vj)= £ £ E Z E Z rap 

in i ——oo m 2 — — co m 3 — — co p\—0 p^— f3= 



7=1 I a ) j=A a j 



x (21) 



xA 



m l m 2 m 3 p 1 p 2 p 3 rm 1 m 2 m 3 p 1 p 2 p 3 i 



where coefficients A mim2m3PiP2P3 are given by (20). 



The case of 



G(xj, vj) = 8( Xl - ft) 5(x 2 - &) 5(x 3 - ft) *( Vl - m) <S(v 2 - « 2 ) <5(v 3 - « 3 ) (22) 



deserves special treatment. Here ft is initial point and w ; is initial velocity. 
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In this case coefficients A 



m 1 m 2 m 3 p 1 p 2 p i 



are equal to 



1 



a \ 



2 1 



x 



xexp 



mi m 2 m }Pl p 2 p 3 2 P> + P2+Pi Pl \p 2 \p 3 \ \2nk) a x a 2 a 3 

Vmim 2 m 3 p i p 2 p 3 (£ji u j)- 



2k {cc^J 



V 



a 2 



J 



(23) 



This leads to following delta function expansion 



d{x x - |0 J(x 2 - £,) <J(x 3 - £,) J(V! - mO <J(v 2 - w 2 ) S(v 3 - u 3 ) = 



a \i 1 ±2? i2? 

~ 7 E E E ex P 

L7ZK ) Cl\Cl2 a i nii=-ca m 2 =-ca m 3 =-oo 



a UrtkX 
2k \ ~cP) 



( 



MJ \ a 2J l fl 3 
V J 



x E E E 

^ 1= p 2 =0 p 3 =0 



<Pm l m 2 m 3 p l p 2 p ) ( x ' ji v j) Ym ] m 2 m i pip 2 p 2 (b p u j) 
2Pi + P2+Plp l \p 2 \p 3 \ 



(24) 



After substitution to (24) expressions (16) and (17), we have 



S(xi - 4i) S(x 2 - | 2 ) <?(*3 - £3) ^(vi - Ui) S(v 2 - u 2 ) S(v 3 - u 3 ) = (25) 



( a \l 1 t2?t£?t2?SSSi^ f (X;-£;)^ f (Vj-UjY 

\^Z^/C y ^^2^3 m,=- 00 m 2 =-oo m 3 =-oo p 1= p 2 =0 p 3 =0 y=l V a jJ\ aa j 

Jzle 2k a2a i ) ( ra~( ATtim ,k Y\ ( rcTf \yiim ,k 



x 



In this way we obtain finally expansion of fundamental solution of Fokker - Planck equation: 

G a (f,x { , v,-, w,-)= (26) 



frl— Z Z Z Z Z £n^W 2 *»'— W 2 *^— V 



Anim :k 
2* I J cr 2 a 



g e 2 H / r^-/ A7tim;k\\ ( rcTf 4mm ,-Jfc" 



We can simplify (26) by explicit summation of the sum on pj. As an instrument for this task we use 
following identity (see [1]). 
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--(x 2 +y 2 ) 

H-OO g 2 

«=0 z " ! 



5" ff„(jc) ff B (y) = 



exp 



2 2/ \2 

x - y (x-ys) 



(27) 



This identity is true when \s\ < 1. 

This gives after long calculations the final form of fundamental solution of initial value problem with 
simple cyclic boundary conditions (4): 



G a (t, x h I,-, v h Ui) = 



1 



a 



3 

2M 



+oo +oo +oo 7=3 

x E E E Ilexp 

W]--oo m 2 --oo W3--00 7=1 



f 2;rra ; 

-fe — '- 

\ acij 



a 1 a 2 a 3 {2ak(l - e~ 2c ") J ^ 
2 1-e" 



nexp 



g (Vj - e" g '« ; ) 2 
2T (1 - e- 2 »0 



Cf 1 + e _£M 



exp 



2;rim f 1 



a 



< (28) 

2 (v ; - e'^r 
a (l+e-»0 



4. Fundamental solution for the unbounded space 

Direct periods a j to infinity and get from (28) fundamental solution for the unbounded space. The 
Fourier series sum is transformed to Fourier integral. Formally, we perform substitutions coj = Inntjla j and 
da>j ~ 2nla j. 



GJ}, x i ,4 i ,v i ,u i ) = 



1 



(In? \ lnk(\-e- 2at ) J %i 



3 

2j=3 



nexp 



+00 +OO +OO 

J 1 .I'M 

—OO —OO —OO 



exp 



[ a ){ a 1+ e~ at 



exp 



a_ (vj - e- at Uj) 2 
2k (l-e- 2at ) 

2 (Vj-e^'u,) 

io)j\ ( Xj - 4j) + - (vj - uj) J 

J[ J J a J J a (1 + e °") 



(29) 



d(o x d(x>2d(OT > . 



We perform calculations (29) using Gaussian integral 



1 +<x 

— f 
2k J 



l 



2AnD 



p 

: e 4D . 



(30) 



For our task D and pj are equal to 



1 2 (v,--«r a 'Kv) 

<v; (v, »,) / / 

J J a J J a (1 + e at ) 



(31) 
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Use another variable D*: 



D* = 



at (l-e- 2a <)-2(l-e- a ') z 
2a~ 4 



(33) 



This variable was introduced in our former work [5] as determinant of quadratic form matrix. We have: 



D = D* 



2ka 



(l-e- 2 «0 ' 



(34) 



Perform transformation and get another form of fundamental solution of Cauchy problem for the 
unbounded space: 



Gooit, X h ft, V,-, Kj) : 



2xk{\ - e- 2a >) J j 



3 

fiexp 



a (Vj-e- at Uj) 2 
2k (l-e- 2at ) 



(35) 



x-| — — I llexp* 



8 I 2nkaD* 



7=1 



1 (1-e" 2 "') 
4fc£>* 2a 



(*j -#;) + - (Vj " My) / — j 

J J a 1 1 a (1 +e 



Join two exponents 



GJf, x h ft, v,-, w,) = 



2^(l-e- 2 "') 



(36) 



1 



xlf^^T fiexp]- 



(1-e" 2 "') 
2or 



(*j " £j) + - (V; " Uj) J I + 

J J a J J a (1 + e °") 



at (1 - e~ 2a ') - 2(1 - e- a 'f (vj - e- a ' U j) 



at\2 



-at,, \2 



+ 



(l-e" 2 «0 



or 



G co (f, x h ft, v h Ui) ■■ 



1 J= 2 

— nx 

{ATtk^u:^ %\ 



(37) 



xexp 



4fc£>* 



(1 - e- 2at ) 



2a 



\( 1 Mv-.-e- at ui) 

- \( X -4)+-( V - Uj ) y jJ 

a \ J 1 a 1 1 J (1 + e °") 



a A 



In order to perform final step, we introduce another set of variables: 



i,- = jc,- - (ft + -i(l 
a 



(38) 
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Vj = v,- - u t e . 



These variables were introduced in [5] as integrals of characteristic equations. Subexpression of (37) 
is equal to 



1 Vi 

a a 



(39) 



and the whole (37) is now 



G x {t, x h Vj, «,) ^ 



1 



(4xk^D*) 3 7=i 



7=3 



4fc£>* 



(1 



2a 



7 a Jil + e- 6 ") 



1 -2 

31 v 7 



■(40) 



We see, that exponent's argument is quadratic form of new variables x, v. 



Very simple calculation brings this expression to final form: 



Goo(t, Xj, 4 h v h Uj) ■■ 



7=3 

nex P < 



{Ajrk-ID^ %i 4-kD 



(l-e- 2a >) „ 2 
2a Xj ' 



(41) 



As expected, this result coincide with obtained in [5]. 



We checked our considerations, comparing result with earlier one obtained by other means. 



5. Calculation of macroscopic parameters 



In this section we calculate averages for Eigenfunctions of Fokker Planck operator. We use the 
method, developed in our work [3]. 

First of all we get expression for Fourier transform M of eigenfunction (16) velocities. We study case 
of nonlimited space and therefore rewrite (16) once again, using substitutions cOj = 2nmjlaj and 



do)j ~ 2x1 a j. 



f[ex P (^,-^)ex P (-^ v^(^(v i + ^)} (42) 



The Fourier transform of eigenfunction (16) is 
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F(t, x h q { ) ■ 



\3/2 ;= 



27ta ) 



nJ2 



j= 3 (2k 

ITH)"^— ] exp(io)jXj) exp 



CO; 



a 



(43) 



(43) is simply equation (19) from our work [1]. We only restored constant multipliers, which we pre- 
viously dropped for simplicity sake. We see, that (43) is generalization of static solution from our work [3]. 
We get the static solution from (43), when all coj and tij are zero. 

We find derivatives of F easily 



dF_ 
dqj 



a V a 



co, 



1i ~ 



(44) 



dq) 



-\<lj + 
a \ a 



Il- 



ea, 



F 



k 

— + 
a 



co , 



a 



(45) 



d 2 F 
dqfiqj 



k I CO; 

<?; + 



a 



a 



qr 



a 



CO ; 

I <?, + — 

a \ a 



co , 



a 



(46) 



when i& j. 



We can write (45) and (46) as a single expression, using Kronecker delta 



d 2 F 



k I co , 
Qi + 

a 



CO; 



Vi- 



ce 



CO; 

a a 



1i ~ ■ 



F - 8;, 



— + — 

a f 



F. (47) 



Now let us return to expansion (21) for solution of Cauchy problem. We use expression (43) for 
Fourier transformation of each term in (21) and get in this way expression for Fourier transform of arbitrary 
solution 



+oo +oo +oo 



M(t, 



I p p p CO CO OO 

(2ny J J J Pi=0 P2= o P3= o 



—oo —oo —oo 



3 

E 

7=1 



k 3 

a LPj + — L °>) 

01 7=1 



(48) 



( k f 2 ^ n (2kY'' 2 
xAmcoMp^p,]^ J UH) Pi y— J expO'^Oexp 



2a\ q ' + 



a 



qj ] dco l dco 2 dco 3 . 
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+oo +00 +00 



dM 1 fff" S » 

z~ = £ £ £ ex P 

dq k (lay J J J =0^ =0^=0 

-OO -OO -OO ' ' ' Z * J 



3 £3 
7=1 a 7=1 



XA 



a) l a> 2 a> i p l p 2 p 3 



k r 



2na ) 



k ( co k 
-\q k + — 1 + 
a I or 



a 



j =} (2k Y //2 
x]^[(-j)" J — exp(ia)jXj)exp 

7=1 V a J 



CO, 



2a ^ 



qj ] dco x dco 2 dco^. 



1 +oo +oo +oo 

d z M 1 r r r « ~ ~ 

= 7+^ ; £ £ £ ex P 



-oo — oo — OO 



k 3 



7=1 a 7=1 

V J 



xA 



a >\< 2 a >iP\PlPi 



( k v /2 



--!«* + — 1+ ^ 



a 7 



ft) t 



— + 



ft^ 
a 



7=1 



2/c 



a 



p/2 



exp(j'o ; Xj) exp 



it 

2a 



qj + 



co, 



a 



CO, 



qj ) dco 1 dco 2 dco 3 . 



d 2 M 1 



+oo +oo +oo 



| p p p OO OO OO 

= 7^- T £ £ £ ex P 



— oo —oo —oo 



( \ 

3 3 
a E /> j + ~2 E ®? 

V J 



xA 



k \ 



3/2 



m 1 m 2 a 3 p l p 2 p 3 \ 



2na J 



a>k), Pk 



q k + — + 



9f 



+ — | + 



oh 
a 



7=1 



2/c 



P/2 



exp(;'ft) -x )exp 



2a a 



] dco l dco 2 dco 3 



when ^ /. 



Density is equal to 



+oo -t-oo +oo 



C C r oo oo oo 

p = (2;r) 3 M = £ £ £ exp 

J J J Pl =0p 2 =0p 3 =0 



— oo — oo — oo 



f \ 

3 £ 3 

« E + ~2 E ®' 

7=1 



« 2 7- 



-11- 



xA 



k f 2 tl 



U(-O p - 



2k 



Pj/2 



exp(;'«,x,)exp 



k (co 



2a V a 



co , 



dcoidco 2 dco 3 



Average velocities are 



+oo +oo +oo 



u k 



j i c c c 

= (27tf — M k = — I I I exp 

O J J J „_n „.-A „.-A 



P J J J p,-0 p-,-0 /?i=0 
-oo -oo -oo 11 11 ' - 



3 £ 3 

« E/ ? j + — E®? 

7=1 « ;=i 



xA 



a> l m 2 ajip l p 2 p 3 



k \ 



3/2 



2;ra J 



;'= 3 (2k'\ Pi ' 2 
x IlH)'' j — exp(/o 7 -x 7 -) exp 
7=1 \ a J 



k ( coj 



2a \ a 



co, 



a 



dcoxdco 2 dco^. 



The current of momentum tensor is 



p p r* OO OO OO 

/«=-(2*) 3 M H =- E E E exp 

Pi=0p 2 =0j> 3 =0 



— oo — oo — oo 



Ae 

ar 2 i= i 



« E/>, + ^2X 



xA 



0> l W 2 0) 3 p l p2Pi 



f— T 



k a Pi 

2 

a co. 



j= 3 (2k Y 1 ' 2 
x IlH) p ' — exp(icojX j)exp 
;=1 V « y 



2a I a 



a 



dco l dco 2 dco 3 . 



Stresses are equal to 



a u = pu k u, - J kl = {2k) M kl - 



M k M, 
Mo 



Hydrostatic pressure is 



P = 



1 V 

-> k=\ 



1 - ,3^( M M kM k 
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6. Calculation of macroscopic parameters of fundamental solution 



Coefficients A 



m 1 co 2 m 3 p l p 2 p 3 



for the case of fundamental solution for unbounded space according to (23) 



and substitutions a)j = Ixnijlcij and doij ~ 2x/aj are equal to 



1 



a \ 



2 1 



ohhhpiph* 2Pt + P2 + Pi Pl lp 2 \pi! \2xkJ ai a 2 a 3 
2k 



xexp 



a 



aiia) 2 o) J pip 2 P3 



(57) 



Expression for M is (see (48)) 



+00 +oo +oo 



M(t, 



j p r> p 00 CO OO 

(zxy j j j Pi=0 P2= o P2=0 



— oo — oo — OO 



f \ 

3 3 

« E Pi + -2 E «>? 



(58) 



1 



a y. 

2P' + P2 + Pip l ]p 2 ]p :i ] [2zk j 



exp 



2k i 

— \coi + (o l 2 + ml 



'2k 



p/2 



x| j riH) Pj '| — I exp(;« i x J ) exp 



1 ( a 3 f 2 3 
(2^[^ko) nex P He-v ,,,)exp 



2^' + ^ 



ft), 



Qj-—\ Y^m^p^pS^m, u m )dcoidco 2 dco 3 = 



\ flfc XjXj + ia(l ~ e ~ a,)1 Xjqj + a I at (1 " e ~ 2a>) " 2 (1 " e ~ a,)1 



The last expression we know from [3]. In this expression we use following definitions: 

6 = 2at -(!-£>-"') (3 -e~ at \ (59) 



and i, - as in (38). 

The first expression in (58) looks rather awkward, but in fact each term in it is a product of three fac- 
tors, each dependent only on one of p y 

Therefore we can write (58) as a product of three independent sums 

However, this decomposition not possible in the common case. We could not expect, that we can 
simplify expressions (48), (52-56) in the common case. 

To prove (58) we use following identity (ref. [4]) 



OO 

£//„(x)- = exp(2xz-z 2 ). 

n=0 " ! 



(60) 
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In our case this give 



Timk \ 



(61) 



exp 



2iak ) ( co 
' '" a 2 a 



exp 



e~ 2a, ^-\q-- 
2a \ a 



This provide instrument to calculate series in (58). The result is : 



^ +oo +oo +00 

M (f , xj, §j, qj, uj) = ^-^ J" J" J" f[ exp 



-co —oo -co " 



-k\ — \t 
a 



exp licojixj-tj-^-) x (62) 



a 



xexp 



a (l-cOjk 
2k { a 2 



exp 



2a x h ' a 



exp 



-le | Uj + 



2icok\( co-, 
~T - — 



exp 



^ , k ( co , 
e~ 2a < — — 



dco l dco 2 dco 3 



We collect terms with equal powers of to ; : 



+00 +OO +00 . - 

M( f ,x^,, ?J , M ,)= (2 ^ J j jftexp 



-co —co —co ' 



fa a 4r 1 2A: k 

-co"\ 1 Ye e 

J ' a 1 2k a 4 2a a 2 a 3 2« 3 



K63) 



xexp 



Uj k 2q 



at, 11 j 2ik 



-2a t 



a a 1 - 



2a 2q > a 



xexp 



2at k j_ 



2a 



qj - le u jqj + e — qj 



dcoidco2dco2 



or 



+oo +oo +oo . 

1 r r r i=3 



M(t,Xj,€j,qj,Uj)= { { jn eXp 



— oo —oo —oo 



ft kt 3 k . 2k -Tjyt k 
-co \ — -+e at — - e 



J {a 2 2 a 3 



a J 



2« 3 



X (64) 



xexp 



a,j\i( X j-tj-(l-e-'»)%- k ^(l-e-<»f 
a a 1 



xexp 



qj z (l-e- 2a ')-ie-°"ujqj 



2a 



dcoidco 2 dco 3 . 



We use following identity to calculate integrals in (64) 



JL f 

In J 



e -Do>\icof d(0 = 



1 



2AkD 



:C 4fl. 



(65) 
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where detailed expressions are 



fj = ( Xj e~ a ') U ^) + i k ^(l- e- at ) 2 . 



a cr 



(67) 



The result of our calculations is 



M(t , xj, 4j, (/,. uj) 



1 Jp3 



(2;r) 3 8(;r£>) 3/2 



1 

4L> 



( Xj - 4 j - (l - ^) + ^ (l - e -') 2 



K68) 



xexp 



2a 



^(l- e - 2 » f )-/ e -»'„. ? . 



Return to "hated" variables (recall (38)): 



M(t,Xj,%j,qj,Uj) = 



1 ( ir r _ i 



3 A 3/2 ,=3 



(2;r) 3 ^2^(9 J %\ 



(69) 



xexp 



2k9 



x+ ikq ^(l-e- a 'f 



2a 



We get following expression, which is identical to last expression of (58). 



M(t,Xj,4j,qj,Uj) 



Y[exp\-ie a, Ujqj X 



(70) 



xexp< 



^ Xj + /«(! - e""') 2 ^ + A [ 2«,(1 - e " 2 -) - 4(1 - e -) 2 1 <? 2 



1 f a 3 V' 2 
= exp(-ie~ at u ,o,)exp 



J + '«(1 " e""') 2 + ^\at(l- e- 2at ) - 2 (1 - e""') 2 ] q) 



In the following we spare such long calculations and provide only their results. 
The density is 
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+oo +oo +00 



(*[*[* oo oo oo 

p = (2;r) 3 M = E Z £ exp 

J J J n.-^ n.-A n^-H 



— oo — oo — OO 



pi=o p 2 =e> ft=0 



3 £ 3 

7 I " 7=1 



(71) 



1 ( a \2 



2k 



2 , 2 , 2 
G) l +6) 2 + a 3 



k f /2 7=3 (2k Y j/2 







2or ^ ar J 





co l m 2 ai 3 p l p2P i 



(4 m ,u m )do) l da 2 do} 3 = 



r 3 x 3/2 

f cr \ (la 



i2»wj exp r*27 w 



Average velocities are 
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The current of momentum tensor is 
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APPENDIX 

We obtained expressions for eigenfunctions (16), (17), and expression for fundamental solution (26) 
for the case of cyclic boundary conditions (4). In this appendix we give formulae for the case of 
unbounded space. They are: 
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DISCUSSION 

In this paper we get the closed form expression for eigenfunctions and eigenvalues of 3D Fokker - Planck 
differential operator. This experssion is known from previous works - see [5]. We use the decomposition 
to obtain solution of Cauchy problem for 3D Fokker - Planck equation. Both periodic problem with cyclic 
boundary conditions and the unbounded space case are considered. We calculate macroscopic parameters 
of 3D Fokker - Planck flow also. 
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